Test of the Equivalence Principle Using a Rotating Torsion Balance 
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We used a continuously rotating torsion balance instrument to measure the acceleration differ- 
ence of beryllium and titanium test bodies towards sources at a variety of distances. Our result 
Aajv,Be-Ti ~ (0.6±3.1) X 10~^^ m/s^ improves limits on equivalence-principle violations with ranges 
from 1 m to cx) by an order of magnitude. The Eotvos parameter is riEarth.Be-n ~ (0.3±1.8) x 10"^"^. 
By analyzing our data for accelerations towards the center of the Milky Way we find equal attrac- 
tions of Be and Ti towards galactic dark matter, yielding rjDM,Be-Ti ~ (—4 ±7) x 10~^. Space-fixed 
differential accelerations in any direction are limited to less than 8.8 x 10~^^ m/s^ with 95% confi- 
dence. 

PACS numbers; 04.80.Cc 



The equivalence of gravitational mass and inertial mass 
is assumed as one of the most fundamental principles 
in nature. Practically every theoretical attempt to con- 
nect general relativity to the standard model allows for 
a violation of the equivalence principle [T]. Equivalence 
principle tests are therefore important tests of unification 
scale physics far beyond the reach of traditional particle 
physics experiments. The puzzling discoveries of dark 
matter and dark energy provide strong motivation to ex- 
tend tests of the equivalence principle to the highest pre- 
cision possible. 

Over the past two decades we have conducted labora- 
tory tests of the equivalence principle 0, i, i, H. This 
letter reports our latest and most precise measurement 
using a new, continuously-rotating torsion balance. The 
torsion balance compares the horizontal accelerations of 
test bodies made from two different materials. Acceler- 
ation differences that depend only on the test-body ma- 
terial violate the equivalence principle. We parameter- 
ize such equivalence-principle violating interactions by a 
Yukawa potential, which for two point objects is 

yW=aG(^),(i),™e— /\ (1) 

where the interaction range A — h/{mhc) is given by the 
Compton wavelength of the presumed exchange boson of 
mass rrib, and which couples to the "new charge" q. The 
coupling strength a is measured in units of the New- 
tonian gravitational constant G, and /i represents the 
mass in atomic mass units. The instrument consists of a 
highly sensitive torsion balance that is continuously ro- 
tated by a turntable. The torsion balance is composed 
of a material-composition dipole pendulum suspended by 
a fine tungsten wire. Any material-dependent accelera- 
tion on the composition dipole pendulum to objects not 
on the turntable produces a periodic twist of the torsion 
fiber. The twist angle, Q, was recorded using a coro- 
tating autocoUimator. The phase and frequency of the 
sinusoidal variation of allowed us to search for differ- 
ential accelerations towards terrestrial and astronomical 
sources. 



Figure [T] shows a schematic drawing of the apparatus 
and the 70.3 g pendulum. The pendulum body is a thin 
aluminum shell with fourfold azimuthal symmetry and 
up down reflection symmetry. It carries four beryllium 
and four titanium test masses in a horizontal dipole con- 
figuration. These two materials were chosen primarily 
to maximize the difference in baryon number {B/fi is 
0.99868 for Be and 1.001077 for Ti), and secondly for 
experimental reasons, such as densities, magnetic prop- 
erties and machinability. The Ti test bodies are hollow 
to match the external shape and mass of the 4.84 g Be 
test bodies to within 50 /ig. The test body shape allows 
us to reproducibly interchange the test bodies, to mini- 
mize alignment errors, and to equalize their gravitational 
interaction. The optical beam of the autocoUimator is re- 
fiected from one of four mirrors located at the pendulum's 
midplane. The entire pendulum and all surfaces near the 
pendulum are plated with « 300 nm of gold. 

The 1.07 m long torsion fiber is a 20 /im thick tungsten 
wire with a torsion constant of k = 2.34 x 10^^ Nm/rad 
yielding a free torsional oscillation period of Tq — 798 s 
and a quality factor Q = 5000 ± 200. The top of 
the torsion fiber is attached to an eddy-current swing 
and vibration damper and tilt isolator. This assembly 
can be rotated with respect to the vacuum chamber by 
a small rotation stage. The vacuum is maintained at 
« 10~^ Pa by an ion pump. The pendulum is surrounded 
by three layers of /i-metal shielding on the rotating plat- 
form and one non-rotating /i-metal shield. To reduce the 
sensitivity to temperature fluctuations and gradients a 
thermally-isolated copper tube inside the vacuum sur- 
rounds the fiber, the aluminum- vacuum chamber is sur- 
rounded by several layers of rotating and non-rotating 
thermal shields and insulation. The entire instrument is 
in a Styrofoam enclosure and the apparatus is located in 
a temperature-stabilized underground room. 

The turntable, below which the torsion balance is at- 
tached, is made from a custom-built aluminum air bear- 
ing 6]. The turntable rotation rate is controlled us- 
ing feedback to an optical angle encoder with two read 
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FIG. 1: Cross section of the apparatus (upper part). The 
entire torsion balance is suspended below a continuously ro- 
tating turntable. Gravity gradient compensator masses were 
placed around the pendulum to reduce coupling to ambient 
gravitational gradients. The pendulum (lower part) carries 
four Ti and four Be masses in a composition dipole. 



heads. The loop is reahzed by a digital signal processor 
controlling an eddy-current drive. The constancy of the 
rotation rate was limited by the angle encoder's linear- 
ity. Angle encoder imperfections were mapped by oper- 
ating the torsion balance at rotation periods Tt <C Tq 
and the corrections were included in the feedback algo- 
rithm. The turntable rotation frequency was normally set 
to LUt/{2TT) ^ 2/3 = 0.835 mHz, where the signal- 
to-noise ratio was found to be optimal. The pendulum's 
angular position and the signals of 31 sensors for temper- 



ature, tilt, etc. were recorded every 2.76 s. 

Gravitational forces between our pendulum and local 
gravitational gradients can occur at the signal frequency. 
To minimize coupling to ambient gravitational field gra- 
dients, Qii, the pendulum was highly symmetrical, with 
the nominal mass moments, vanishing for I < 7. 
In addition, the pendulum's q2o, qso and 540 moments 
were designed to be zero to avoid gravitational coupling 
due to a small misalignments of the pendulum The 
ambient Q21 and Q31 fields at the pendulum's position 
were compensated with 888 kg of lead and 8.8 kg of alu- 
minum, respectively. Lead bricks placed about 1 m from 
the pendulum were used to null the remnant (521-field. 
The gradient fields were measured with a special gra- 
diometer pendulum that was configured to have either a 
large q2i or q^i moment. We found that the (521-field 
varied by as much as ±1% during a year, which we at- 
tribute to changes of the water table. Once the equiva- 
lence principle pendulum was installed, its residual (721- 
and gsi-moments were measured by rotating the Q21- and 
Qsi-compensators by 180°, doubhng the uncompensated 
fields. The pendulum's unwanted mass moments were 
then minimized by adjusting four screws on the pendu- 
lum body until only a small gravitational coupling re- 
mained (see Table HI, which was later subtracted from 
the data. 

The turntable must rotate about local vertical since a 
tilt of the attachment point of the torsion fiber causes a 
small apparent rotation of the pendulum. The tilt of the 
turntable was continuously measured with level sensors 
on the rotating platform. The once-per-revolution com- 
ponent of the tilt was minimized by a feedback loop that 
changed the length of two of the turntable's support legs 
by controlling their temperature with Peltier elements. 
This system nulled the periodic tilt of the level sensor 
to within ±3 nrad. At the pendulum body, 1.7 m be- 
low the feedback sensor, local vertical was different by 
55 nrad and a correction to the data was required. We 
inferred this tilt by using a second tilt sensor 0.2 m below 
the pendulum. We also found the change in local verti- 
cal consistent with our local mass integration. The data 
were corrected using the tilt at the pendulum and a tilt 
matrix characterizing the tilt sensitivity. The tilt matrix 
was measured for all four pendulum mirrors by deliber- 
ately tilting the turntable rotation axis. The magnitude 
of the tilt matrix ranged from « 1% to « 7% depending 
on the mirror. 

The largest entry in our error budget arises from sys- 
tematic and statistical uncertainties associated with tem- 
perature gradients and fluctuations. We assessed the sen- 
sitivity to temperature gradients by placing temperature- 
controlled panels on opposite sides of the apparatus to 
produce exaggerated thermal gradients. The quoted un- 
certainty is limited by sensor noise in resolving the small 
temperature gradients that occurred during equivalence 
principle data taking. By applying a magnetic field and 
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field gradients using coils, as well as strong permanent 
magnets, we found the magnetic coupling to be small 
(Table m. 
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FIG. 2: Sliown are measured differential accelerations to- 
wards North (top) and West. After tfie first four data runs 
tfie Be and Ti test bodies were interdianged on the pendulum 
frame. A violation of the equivalence principle would appear 
as a difference in the means (lines) of the two data sets. The 
offset acceleration is due to systematic effects that follow the 
pendulum frame but not the composition dipole. The data 
have been corrected for tilt and gravity gradients, but only 
the statistical uncertainties are shown. 



with the torsion fiber, the pendulum body, or the mag- 
netic damper. An equivalence principle violation would 
cause a difference in the measured acceleration of the two 
sets (Fig. [2]). 

Our data analysis began by using a digital notch fil- 
ter to remove the pendulum's free torsional oscillation. 
The data were then divided into segments containing 2 
turntable revolutions (w 2400 s). Nine harmonics of the 
turntable rotation frequency, an offset, a linear and a 
quadratic drift were fitted to the data in each segment. 
About 7% of the segments were excluded from the analy- 
sis due to spikes in the ion pump current or due to abrupt 
changes of the turntable axis from local vertical. The co- 
efficients of the sin (tott) and cos (cott) terms were used to 
extract the signals of interest and their scatter was used 
to determine the statistical uncertainty. The pendulum 
data were corrected in amplitude and phase for the pen- 
dulum's dynamic response to external torques and for 
electronic attenuation and time delays. After correcting 
for systematic effects (Table U), the data were reduced to 
material-dependent differential acceleration components 
along North and West directions. Figure [2] shows these 
components for six data runs. We observed a small unex- 
plained offset in the East- West direction that followed the 
suspension system (pendulum body, torsion fiber, swing 
damper) but the signals tracking the composition dipole 
were unresolved in all directions. 

The composition dependent accelerations are 

aAr(Be) - aA,(Ti) = (0.6 ± 3.1) x 10"^^ m/s^ and 
aw{Be) - aw{Ti) = (-2.5 ± 3.5) x 10"^'^ m/s^ 



TABLE I: The raw differential accelerations between Be and 
Ti towards North (N)and West (W) are shown in line 1. Lines 
2 to 5 list corrections that were applied and the bottom line 
gives our corrected results. Uncertainties are la. 
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We collected 75 days of equivalence-principle data in 
two sets. The orientation of the pendulum with respect 
to the vacuum chamber was changed by 180° once a day, 
mostly to cancel the effect of small turntable rotation 
rate imperfections. Approximately biweekly we changed 
the orientation of the pendulum by 90° and took data 
using opposite mirrors. Between the two sets, we inter- 
changed the test bodies with those on the other side of 
the pendulum, reversing the composition dipole on the 
pendulum frame to eliminate spurious signals associated 



with titanium being more attracted to the South and to 
the West. Our best limit on the classical equivalence 
principle parameter [3| t] is 



77(Be - Ti) = 



(0.3 ± 1.8) X 10 



-13 



(2) 



Figure [3] shows the limits on the strength, a, of a new 
interaction (Eq. 1) as a function of range A. To estab- 
lish these limits we used the mass density and charge 
content of the environment surrounding the torsion bal- 
ance to create a source model. For A = 1 — 100 m the 
source is dominated by a hill sloping towards the East. 
For A < 10 km the local topography and bedrock become 
significant. At ranges between 10 km and 1000 km, pre- 
liminary results using large scale density and composi- 
tion models indicate that the limit on a is better than 
the dashed line shown in Fig.[3l A detailed description of 
the model and limits will be included in a future publi- 
cation. We used an elliptical layered Earth model[3,[l, Q 
for A > 1000 km. For this range the source mass is lo- 
cated towards the North. 

Equivalence-principle violating interactions associated 
with an astronomical source are additionally modulated 
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by a solar or sidereal frequency: 

Aajy— cos0( — Aacos (0 — (jjo) — Aasin (0 — 0o)) +on, 
Aaw— cos9(^Aa sin {(p — ^q) — Aa cos (0 — (j>o)) +ow, 

with and </> being the altitude and azimuth of the astro- 
nomical source, Aa the differential acceleration towards 
the source and Aa its quadrature component; on and 
ow are possible instrument offsets. Figure H] shows the 
averaged AaAr and Aa^y versus sidereal time. A simul- 
taneous fit of AaAT and Aaw towards the galactic center 
yields 
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a{Be) - a{Ti) 
h{Be) - h{Ti) 



Aa 
Ad 



(-2.1 ±3.1) X 10"^^ m/s^ , 
(2.7 ±3.1) X 10"^^ m/s^. 



Since only about a quarter of the total acceleration of the 
solar system towards the center of our galaxy is caused 
by galactic dark matter [l3|, we find rj]jM.Be-Ti — (— 4± 
7) X 10"^ 

With 95% confidence we constrain space-fixed differ- 
ential accelerations in any direction to be smaller than 
Aa = 8.8 X 10-1^ m/s^. 
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FIG. 3: New upper limits on Yukawa interactions coupled to 
baryon number with 95% confidence. The uncertainties in 
the source integration is not included in this plot. The num- 
bers indicate references. The shaded region is experimentally 
excluded. Preliminary models for 10 km < A < 1000 km 
indicate that the Umit on a is smaller than the dashed line. 

We have substantially improved the limits on the 
strength of an equivalence-principle violating, long-range 
interaction. We are currently broadening our search by 
using other test-body materials and improving the sensi- 
tivity of our torsion balance. 
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FIG. 4: The averaged differential acceleration of Be and 
Ti towards North and West as a function of sidereal time. 
The dashed line represents a hypothetical signal of 20 x 
10~^^ m/s^. The solid line is the best fit toward the galactic 
center (Aa = (-2.1 ± 3.1) x 10"^^ m/s^). 



Schmidt, Tom Butler and Chris Spitzer have contributed 
to the apparatus development. 
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